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ABSTRACT 

Coronal mass ejections (CMEs) are powered by magnetic energy stored in non-potential (current-carrying) coronal 
magnetic fields, with the pre-CME field in balance between outward magnetic pressure of the proto-ejecta and 
inward magnetic tension from overlying fields that confine the proto-ejecta. In studies of global potential (current- 
free) models of coronal magnetic fields — Potential Field Source Surface (PFSS) models — it has been reported 
that model field strengths above flare sites tend to be weaker in when CMEs occur than when eruptions fail to occur. 
This suggests that potential field models might be useful to quantify magnetic confinement. One straightforward 
implication of this idea is that a decrease in model field strength overlying a possible eruption site should correspond 
to diminished confinement, implying an eruption is more likely. We have searched for such an effect by post facto 
investigation of the time evolution of model field strengths above a sample of 10 eruption sites. To check if the 
strengths of overlying fields were relevant only in relatively slow CMEs, we included both slow and fast CMEs in 
our sample. In most events we study, we find no statistically significant evolution in either: (i) the rate of magnetic 
field decay with height; (ii) the strength of overlying magnetic fields near 50 Mm; (iii) or the ratio of fluxes at low 
and high altitudes (below 1.1R@, and between 1.1— I.SRq, respectively). We did observe a tendency for overlying 
field strengths and overlying flux to increase slightly, and their rates of decay with height to become slightly more 
gradual, consistent with increased confinement. The fact that CMEs occur regardless of whether the parameters we 
use to quantify confinement are increasing or decreasing suggests that either: (i) the parameters that we derive from 
PFSS models do not accurately characterize the actual large-scale field in CME source regions; or (ii) systematic 
evolution in the large-scale magnetic environment of CME source regions is not, by itself, a necessary condition 
for CMEs to occur; or both. 
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1. INTRODUCTION 

In a coronal mass ejection (CME), the Sun launches a hot (~ 1 
MK) mass (~ 10 15 g) of magnetized plasma at high speed (from 
a few hundred km sec -1 at the slow end to 2000 km sec" 1 or 
more at the fast end) into interplanetary space. A solar flare 
is a transient burst of energy released across a wide range of 
the electromagnetic spectrum, from radio to X-ray (and in some 
cases, gamma-ray) wavelengths, over a wide range of energies, 
from ~ 10 - 10 33 ergs for X-class flares at the high end to 
barely detectable X-ray enhancements in weak A-class flares. 
CMEs are often associated with flares, but the relationship be- 
tween the two is not one-to-one: some flares produce no CMEs; 
some (mostly slow) CMEs occur with barely any detectable 
flare emission; but most fast CMEs are associated with flares, 
and almost all large (X- and M-clas s) flares produce CMEs 
( Andr ewilOOl iWang & Zhandl2007h . 

CMEs are observed to originate in the low corona. Considera- 
tion of available sources of energy there imply that both flares 
and CMEs mu st be powered by energy stored in the coronal 
magnetic field (Forbes 2000): coronal gas densities are so low 
that thermal energy, gravitational energy, and kinetic energy 
densities are all too small to supply the energy required to drive 



CMEs. 

The basic physical processes that produce CMEs, however, are 
not understood. Beyond our scientific interest in CMEs, those 
directed toward Earth are the primary drivers of severe space 
weather events. While flares also occasionally produce effects 
at Earth independent of CMEs — for instance, ionization of 
the Earth's upper atmosphere from the enhanced X-ray flux can 
interfere with radio commu nications — the biggest terrestrial 
effects arise from CMEs (Goslin3 ll993h . Consequently, the 
practical capability to predict CMEs would be quite useful. Be- 
cause the physical processes at work in CMEs are not under- 
stood, however, we are still at the stage of seeking empirical 
correlations between solar observations and the occurrence of 
CMEs. 

The CME process can be understood in terms of the interplay 
between magnetic pressure and tension forces in the corona. 
(The same reasoning used to argue that CMEs must be pow- 
ered by magnetic energy also implies that, in the corona, non- 
magnetic forces — such as gas pressure gradients and gravity — 
are not the primary drivers of CME dynamics.) The pre-CME 
field is thought to be in balance between outward magnetic pres- 
sure and inward confinement by magnetic tension. 
Essentially, CME initiation requires an increase in outward 



magnetic pressure, or a decrease in confinem ent, or both. In the 
tether-cutting model of Moo re et al.1 (120011) . for instance, field 
lines that arch over the proto-ejecta can reconnect underneath 
it, converting downward tension into upward tension, which 
both increases the outward force and decreases the inward force; 
therefore the system will expand outward. As more underlying 
field lines reconnect, the outward forces increase. If the increase 
in outward forces is insufficient to overcome the confinement 
by magnetic tension, then no eruption occurs. If, however, the 
outward forces overcome the restraini ng tension, then an erup - 
tion occurs. In the breakout model of lAntiochos et al.1 _1999), 
"breakout reconnection" above the proto-ejecta decreases over- 
lying magnetic tension, causing the proto-ejecta to rise, which 
initiates "flare reconnection" underneath it, leading to CME ini- 
tiation. 

In contrast to these non-ideal models of CME i nitiation, ideal 
MHD instabilities might also produce CMEs. Kliem & Ttirok 
(2006) suggested the torus instability could explain the CME 
process: if the outward force from magnetic pressure in the 
proto-ejecta decays with radial distance from the photosphere 
more gradually than the inward magnetic tension of the overly- 
ing field that confines the proto-ejecta, then any outward pertur- 
bation of the proto-ejecta will be unstable, and an ejection will 
occur. In a given configuration, variations in either the field of 
the proto-ejecta or the confining field (or both) could enable the 
torus instability to occur. 

One key property of flares and CMEs is that they occur along 
polarity i nversion lines (PILs) of the photospheric magnetic 
field (e.g.. iFalconer et aTll2006l) . Observers report that the ori- 
entations of H-a, EUV (extreme ultraviolet), and SXR (soft 
X-ray) structures low in the solar atmosphere above PILs typi- 
cally deviate from directions expected from potential field mod- 
els, implying electric currents are present; but higher in the 
corona, EUV and SXR loop s typically appear nearly potential 
dMartin & McAllisteri[T996T) . 

These observations of potential-like field configurations far 
above PILs suggest potential field models could be studied to 
understand the role o f magnetic field confinement in CMEs. 
IWang & Z hang (2007) studied potential field models associated 
with eight X-class flares, four with and four without CMEs, 
computing the ratio of flux at low altitude in the corona, be- 
low I.IRq, to the flux from I.IRq - 1.5Rq for each case. They 
found flares with less confinement relative to magnetic pres- 
sure (a higher ratio o f low flux t o high flux) were more likely 
to erupt. In addition, iLiul d2008) analyzed several failed erup- 
tions and CMEs, fitting the decrease with radius of horizontal 
field strength in global potential models, and found that field 
strengths decreased faster with radius in successful eruptions. 
He also found average field strengths at a height of 42 Mm were 
higher in cases with failed eruptions. 

Here, we also investigate the magnetic structure above eruptive 
PILs in potential magnetic field models, to look for patterns in 
the evolution of the large-scale model fields around the times of 
CMEs. In the next section, we discuss how we chose our events, 
then we discuss our methods to analyze the data in |3] look for 
evidence of systematic magnetic field evolution within the data 
in §|U an d finally discuss our results in §5\ 

2. DATA 

For this research, we selected 10 CMEs over a 13 year span, 
from 1997 to 2010. Our goal was to compile an event sam- 
ple commensur a te in l ength with those used by Liu (2008) and 
Wang & Zhang] (120071) . To select the events we analyzed, we 
first manually identified source active regions (ARs) for 100 



CMEs that were within a box measuring 60° by 60° at disk 
center. CME dates and times were taken from SOHO/LASCO 
CME catalog The RHESSI B and NGDC X-ray flare listfQ 
combined with University of Hawai'i's Institute for Astronomy 
active region map£| (compiled from Solar Region Summaries 
released by NOAA's Space Environment Center), were used to 
identify CME source active regions. The requirement that CME 
source regions be unambigously identifiable in data from the 
CME and flare catalogs was a restrictive criterion. Starting from 
the CDAW CME catalog, we recorded each usable CME event 
until we found 100. This compilation was not meant to be ex- 
haustive; we sought, instead, only a sample of CMEs for which 
source active regions had been identified. 
We then distilled this list down to 10 by selecting the most iso- 
lated, compact, and near-disk-center events. We chose ARs 
which were spatially isolated as well as compact; these con- 
straints reduce the likelihood that our magnetic field models 
will be strongly affected by neighboring active regions. We also 
chose a sample that includes both fast and slow CMEs to see if 
the magnetic environment of CME source regions plays a strong 
role in either fast or slow populations. We defined a fast CME 
to be above 800 km sec -1 ; because most CMEs are slower than 
this, we have more slow CMEs in our sample. The 10 ARs 
which best fit these criteria were used for modeling in this re- 
search. Basic properties of our events are presented in Table 
Q] Our events consist of generally weaker flares than those of 
Wang & Zhang (all X-class) and Liu (several X- and M- classes, 
in addition to a few C-class flares). We note that the May 1997 
event has been stu died previo usly by several researchers, e.g. 
iLi et al.1 d2010l) and lLlul d2008l) . 

3. METHODS 

We build upon previous work by looking for any sys t ematic 
time evolution in parame ters like those studied by ILiul ([2008 ) 
and lWang & Zhangl (120071) prior to CMEs. These authors com- 
pared CMEs and failed eruptions in different active regions, 
and found that in potential field models of active regions that 
produced CMEs (cf., failed eruptions), either overlying field 
strengths decreased more rapidly with height or there was less 
total flux overlying the PIL. These authors suggest both that the 
confinement by magnetic fields above active regions is a key 
factor in determining whether an eruption can occur, and that 
potential field models can usefully quantify this confinement. 
Neither paper discussed any prediction about relationships be- 
tween time evolution in overlying fields and the occurrence of 
a CME within an active region. We conjecture, however, that 
if overlying fields do in fact play a major role in confining pre- 
CME fields, then it is reasonable to expect overlying fields to 
weaken around the times of CMEs. The occurrence of CMEs 
despite lack of weakening in overlying fields in potential mod- 
els would suggest that either (i) confinement by overlying fields 
is not a key factor in CMEs; or (ii) potential field mode ls do 
not accurately quantify confinement; or both. Liu (2008) also 
mentions that the behavior of transverse magnetic field strength 
at low heights may also play a role in CME onset. We test our 
conjecture in our sample of CMEs, to see if these parameters 
of the CME magnetic environment change systematically from 
about 48 hours before to 24 hours after each CME. 

1 http://cdaw.gsfc.nasa.gov/CME_list/UNIVERSAL/text_ver/univ_all.txt 

2 http://hesperia.gsfc.nasa.gov/hessidata/dbase/hessi_flare_list.txt 

3 ftp://ftp.ngdc.noaa.gov/STP/space-weather/solar-data/solar-features/solar- 
fl ares/x-rays/goes/ 

4 http://kopiko.ifa.hawaii.edu/ARMaps/Archive/ 
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Table 1 

CME and Active Region Table 



NOAA 


Lat./Long. 


CME Date 


Flare Time, 


CME Time 


CME Speed 


GOES Flare 


B t at t = b 


ARNo. 


(cleg) 


yyyy/mm/dd 


t = (UT) 


(UT) 


(km s" 1 ) 


Class 


(G) 


8038 


N21W09 


1997/05/12 


05:26 


05:30 


464 


C1.3 


9 


10311 


S12E04 


2003/03/13 


02:13 


02:54 


1021 


C1.3 


to 


10696 


N09E05 


2004/1 1/06 


00:57 


01:31 


818 


M5.9 


30 


10715 


N05E06 


2005/01/03 


04:22 


05:06 


363 


C3.8 


7 


10775 


N06W07 


2005/06/12 


02:36 


02:36 


590 


C3.5 


2 


10822 


S07E02 


2005/11/19 


23:38 


00:30 


231 


C4.1 


35 


10956 


N03E00 


2007/05/19 a 


13:02 


13:24 


958 


B9.5 


10 


10960 


S07W02 


2007/06/07 


17:20 


17:30 


263 


Cl.l 


40 


10977 


S05W07 


2007/12/07 


04:41 


05:30 


209 


B1.4 


2 


10978 


S09W10 


2007/12/12 


05:32 


06:30 


305 


B2.4 


40 



d While the CME appeared in LASCO on November 19, the flare occurred the previous day. 
k The mean B, at r = 46.5 Mm, averaged over the photospheric active-region mask; see !4.2l 

We investigate model fields up to 24 hours after the CME event 
for two reasons. First, we wanted to know if the evolution 
of the model field around the time of the CME was part of a 
longer-term trend. If a CME occurs in the middle of an episode 
of persistent flux emergence, but there is no particular tem- 
poral structure in the evolution of the large-scale model field 
around the time of the CME, then we have learned something: 
while flux emergence might destabilize the coronal field to pro- 
duce CMEs, this property of the model potential field cannot be 
used to detect such destabilization. Second, the additional post- 
CME model coronal fields are based upon independent mea- 
surements, and can therefore reveal whether changes in model 
field parameters near the CME are temporary fluctuations, and 
might therefore be ascribed to noise. 

To model the magnetic environment of CMEs, we use global 
Potential Field Source Surface (PFSS) models pro duced with 
the codes developed by Schrijver & De Rosa (2 003b from mag - 
netograms observed by SOHO/MDW (IScherrer et all U995). 
This model's photospheric fields are hybrids of recent and past 
observations, combined with a flux transport model. Recent 
96-minute-cadence magnetograms are used for the visible disk, 
and the flux transport model includes diffusion and a statistical 
treatment of small-scale flux emergence for the far-side bound- 
ary condition. Except where explicitly stated otherwise, we use 
PFSS models computed with spherical harmonic coefficients up 
to r? max =192. The models were computed 6 hours apart, max- 
imizing time resolution while minimizing information redun- 
dancy. The photospheric Alfven speed — of order 1 km s" 1 — 
governs the timescale of photospheric magnetic evolution. In 
contrast to the spatial resolution of the magnetograms assimi- 
lated into the models (about two 1 .4 Mm pixels at disk center), 
the spatial resolution of these PFSS models is about two he- 
liocentric degrees, or ~ 24 Mm near disk center, implying a 
photospheric Alfven crossing time for model pixels of about 6 
hours. (In the corona, the Alfven speed is of order 1000 km s , 
implying dramatically faster evolutionary timescales.) Since the 
potential field models were only constructed every 6 hours, we 
matched the CME event times to the time of the closest model 
field, which we define to be t = to = 0. 

Since we are interested in how our results compare to iLiul 
(2008), we followed his approach and fitted transverse field 
strengths over a range of radii similar to his, at heights of {46.5, 
59.2, 72.3, 85.8, 99.9, 114.5} Mm, assuming the transverse 
field is described by a power law in this range. Transverse field 

5 http://soi.stanford.edu/magnetic/index5.html 



strength is defined as: 

B, = yfBl + B 2 g , (1) 

where and Bg are the orthogonal components of the mag- 
netic field parallel to the Sun's surface. 

We are interested in the time evolution of the magnetic field 
decay rate, 7, as a function of height, henceforth assumed to 
follow a power law in radius with exponent 7, where 

B t oc f< , (2) 

and r is the distance above the solar surface. A more negative 
value corresponds to a magnetic field which is decreasing in 
strength more steeply with radius. In order to compute the best 
estimate of 7, we tried two approaches: averaging the magnetic 
field at each height within an AR then fitting a 7 to the averages 
(Method 1); and calculating a 7 for each pixel within an AR 
and then averaging the 7's (Method 2). Furthermore, we apply 
Methods 1 and 2 for two separate areas: the entire AR itself, 
as well as just the PIL. We analyzed the model fields near PILs 
separately because PILs are source regions for CM Es. 
We use a linear, least-squares fit dPress et al.lll992l) to the log- 
log values of magnetic field strength over radial height to find 
power-law exponents. When fitting the logarithms, we set the 
origin of the independent variable, log(r /Vo), to be between the 
third and fourth radial values. This enables us to plot lines de- 
picting the fitted 7's with errors more clearly without reference 
to variations in the fitted y-intercept, which we did not investi- 
gate. 

For Method 1, uncertainties on 7 are derived from the standard 
error in the estimate of the average magnetic field at each height 
through the fitting procedure. To estimate the uncertainty in the 
average transverse field at a given radius r,, we first compute 
the standard deviation erg, of B, at r, over pixels corresponding 
to the photospheric mask (AR or PIL), and then divide cr» ,■ by 
^//Vpix, where ^N p - lx is the number of pixels in the mask. This 
treats B, from each pixel in the mask as an independent estimate 
of the mean field strength at r, . An array of these uncertainties, 
one per each fitted radial point, is passed to the least-squares 
fitting procedure, where it is used in weighting the points to be 
fit. The procedure computes the square of the uncertainty in 
fitted slope 7 as one over the variance of {ri/(aBtj/ \/N P i X )}, 
multiplied by the sum over i of {1 /(ctbij/ \/Npix)}- For Method 
2, the uncertainty in 7 at a given time is the standard error in 
the estimate using the fits to 7 from each pixel: the standard 
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deviation of the set {7/}, divided by \/Npix, where jj is the 
linear fit to B,(r) in the y'-th pixel, and N p i x is the number of 
pixels in the appropriate mask. 




175 180 185 190 195 200 

Degrees Longitude 



Figure 1. Pixels used in the June 07, 2007 AR. The red and green outlined re- 
gions represent the AR and PIL mask, respectively. The background grayscale 
shows the radial magnetic field (white positive, black negative) at the photo- 
sphere in the PFSS model. The AR mask used in finding 7 is defined to be 
the pixels whose absolute radial field strength is above 25 Mx cm - -. The PIL 
mask is computed by dilating the masks of positive and negative polarities with 
absolute radial field strengths greater than 50 Mx cm -2 by 1 pixel and finding 
where they overlap. The grayscale saturation level of the image is set to ± 50 
Mx cm -2 . 

Figure 03 shows how we mask the AR and select pixels from 
which to compute 7. The background grayscale shows the 
radial magnetic field from the PFSS model at r = 0, with 
white representing positive fields and black representing neg- 
ative fields. We computed the AR mask by identifying pix- 
els with an absolute radial magnetic field strength above 25 
Mx cm" 2 . We iden tified PIL masks fo llowing the method of 
lSchriiverl(l2007l) and lWelsch & Lil (12008b : masks of pixels from 
each polarity stronger than a threshold are dilated by a pixel and 
multiplied together; regions of overlapping pixels are thereby 
identified as the PIL mask. A higher threshold, 50 Mx cm" 2 , is 
used in finding the PIL mask; starting with too low a threshold 
leaves too many pixels in strong field regions of either polarity. 
Similar magnetograms of all ARs in our sample are shown in 
Figure 2. 

When choosing to analyze the entire AR, the spread of magnetic 
field strengths is typically very wide (see panel (a) of Figure |3). 
Method 2 nearly always results in a 7 with a smaller magnitude 
than Method 1 . In addition to this, it also results in larger uncer- 
tainties in 7. While Method 1 fits closer to the average magnetic 
field, this does not necessarily mean it is a better estimate of the 
decay rate in an AR's overlying magnetic field. Method 2 takes 
into account field data which stray from the average. While the 
individual points have an insignificant impact on the final result 
of Method 1, they appear in Method 2 by shifting the final 7 
down. Each method offers a slightly different perspective on 
field structure. We applied two methods of calculating 7's in an 
attempt to test the robustness of our computation. As shown in 
Figure @ we found that the main factor which determined the 
value of j was not the computation method, but instead which 
part of the AR we looked at. 



Table 2 

Evolution of 7 



Initial Time" 


Final Time" 


A 7A R > b 


A7AR < 


A7pil > 


A7pil < 


-6 





6(1) 


4 


7(1,1) 


3(1) 


-12 





6 


4(1) 


5(1) 


5(1) 


-24 





7(1) 


3 


8(2,1) 


2(1) 


-48 





7(U) 


3 


7(2,1) 


3(1,1) 


-12 


12 


4 


6(1,1) 


3 


7 (2,1) 


-24 


24 


6(1,1) 


4 


6(2,1) 


4(2,1) 


-48 


24 


7(1,1) 


3(2) 


5 (2,1) 


5 (2,1) 



a Measured in hours since closest magnetogram to CME event 

k Most values are not significant at a la confidence level; the number of events for which 
values are significant are listed in parentheses with the first and second entries correspond- 
ing to lcr and 2a significance, respectively. 



As a caveat, due to lack of resolution in the magnetic models, 
for a few events there are very few PIL pixels; in some, there 
are less than 10 pixels in the PIL masks. 

4. RESULTS 

4.1. Evolution in 7 

The purpose of this research is to see if there is any relation 
between evolution of 7 and a CME event. Figure |4] reveals no 
strong pattern in 7's: in most cases, 7 remains relatively un- 
changed, with no strong tendencies (relative to error bars) for 
increases or decreases. We note th at our values for 7 are similar 
to those determined by Liu (2008), who found 7's in the range 
of 1.51 - 2.25 over a range of r adii (42 - 105 Mm) similar to 
ours. (In his Fig. 3, Liu 2008 shows a faster fall-off of field 
strength with radius beyond about 400 Mm.) For June 07 2007, 
we saw a slight increase in 7 (implying a more gradual decay 
of transverse field with radial height), while for June 12, 2005 
we saw a slightly stronger decrease in 7 (steeper decay of field 
strength with radius). (Alone among our events, the Decem- 
ber 07 event exhibits strong back-and-forth fluctuations.) A key 
point to keep in mind is that CMEs can evidently occur for any 
type of evolution in 7. Furthermore, we found no significant 
correlation between j(to) and either CME speed or flare size. 
For the entire set of events considered as a whole, we do not 
see evidence of systematic decreases in 7's around CME onset, 
regardless of whether 7 is computed with Method 1 or 2, or for 
whole ARs or only PILs. 

We note that the 7's computed from PILs are consistently more 
negative than their AR-computed counterparts. One possible 
explanation for this is due to the nature of PILs themselves. 
The PIL lies in the region of the photosphere between the two 
polarities, meaning regions of corona above PILs are farthest 
away from photospheric source regions, resulting in weak po- 
tential fields (which scale as the square of the inverse distance). 
However, this behavior only affects the value of 7, and not its 
evolution over time. 

To investigate whether or not there exists a more subtle pattern 
in the evolution of 7, in Table [2] we show changes in 7 deter- 
mined using Method 1, for several time intervals near each CME 
event. We chose four time interval s prior to CMEs to loo k for 
systematic, pre - CM E changes in 7. lLongcope et al.l (120051) and 
Schr iiver et alj d2005l) both found characteristic time scales of 
coronal energy storage of about 24 hours, suggesting this might 
be a characteristic time relevant to CMEs. Our pre-CME inter- 
vals were chosen to sample near this characteristic time scale. 
In addition, we also computed field changes over one, two, and 
three day intervals that straddled each CME. We find that most 
changes in 7 are not statistically significant at the 1-ct level; 
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Figure 2. Magnetograms of source ARs. The red and green outlines correspond to pixels used for the entire AR and the PIL-region, respectively. Grayscale saturation 
levels are ±50 Mx cm -2 . 
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(c) Fit and Uncertainties in Method 1 



(d) Fit and Uncertainties in Method 2 



Figure 3. Comparison of different Methods 1 and 2 used to calculate 7. In (a), the black dots represent transverse magnetic field strengths at each radius, and the 
average values ( X ) are fitted. In (b), the dashed lines represent the fitted power laws for field strength above each pixel, and the solid line is the average of the fitted 
slopes. We look at values of the fitted power-law exponents (the fitted slopes) of the solid lines over time for each event in our sample. In (c) and (d), we plot fitted 
slopes varied by the estimated uncertainties. In all four panels, the X denote the average field. The boxes in (b) and (d) represent the median of the magnetic field 
values at each radius. While the distribution of field values is wide, the median is close to the average. 



changes that are significant at the 1- and 2-cr level are shown 
in parentheses. We only show values of 7 from Method 1 in 
Table |2]because both methods had similar results for most time 
intervals. 

Table [2] reveals that, while most of the changes in 7 are not 
statistically significant, for A7AR, there is a preponderance of 
positive changes in 7ar: for six out of seven intervals, posi- 
tive A7ar's outnumber negative A7ar's. To test whether this is 
significant or not, we compared our result with the null hypoth- 
esis: that changes in 7ar are random. Since the evolution of 7ar 
is observed to be either positive or negative, the signs of 7ar's 
under the null hypothesis can be modeled with the binomial dis- 
tribution. We first calculated the probability that a positive bias 
in A7AR (more than five events with a positive change) would 
occur in one time interval (0.37); we then calculated the likeli- 
hood that the same bias would occur in at least six out of seven 
intervals (0.01). Consequently, we argue that our results are in- 
consistent with random changes in 7ar (which corresponds to 
the signs of A7 also being random): 7 tends to become more 
positive, implying less rapid decay of B, with height as time 
goes on. There is also a predominance of positive A-fpn's in 
Table|2] but not as great as that among the A7ar's. 
We note that, because the same magnetic model appears in mul- 
tiple time intervals, a few random excursions in 7 at the times 
that we selected to difference 7ar could have introduced a spu- 
rious pattern into our results. We therefore evaluated A7AR for 



all possible time intervals in our dataset. In Figure [5] we show 
positive changes in blue and negative changes in red, for each 
event. The saturation in the color scale is set to ± 2, where 
A7AR for each event is normalized by the standard deviation 
in 7ar for that event. Horizontal or vertical bars in the shad- 
ing would imply magnetic models at one particular time could 
be responsible for the biases seen in Table |2j while some bars 
can be seen, they do not appear to underlie the predominance of 
positive A7AR values in this Table. The smoothness in the shad- 
ing for most events demonstrates that typical evolution in 7ar is 
steady. Looking at all 10 panels, there is a slight predominance 
of blue, consistent with the bias seen in A7AR in Table [2] 
Recall that a positive A7 corresponds to a less-steep decrease 
of magnetic field strength with radius. Decreasing confinement 
would correspond to 7 becoming increasingly negative prior 
to a CME. However, only one event (June 12, 2005) clearly 
showed this behavior. 

4.2. Evolution in B, 

iLiul (120081) also suggests weaker transverse field strengths at 42 
Mm more easily enable CMEs. To test if changes in B, are re- 
lated to CME onset, we also looked at evolution in transverse 
field strength at 46.5 Mm, the height in our models closest to 
42 Mm, which Liu used in this aspect of his modeling. Look- 
ing at evolution of B, plotted for each event in Figure [6] we see 
no typical trend in the evolution of field strengths surrounding 
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Figure 4. Time evolution of 7's. The solid line is Method 2 (AR), the dotted line is Method 1 (AR), the long dashed line is Method 2 (PIL) and the short dashed line 
is Method 1 (PIL). For most events, Methods 1 and 2 give similar answers for AR, and for PIL — but PIL 7's are significantly more negative than AR 7's. A more 
negative 7 corresponds to more rapid decay of the magnetic field with height. Some error bars and data points are slightly offset in time to avoid overlap. The solid 
vertical line corresponds to t = and the dashed vertical line corresponds to the CME time. 
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Figure 5. A7 between all possible pairs of 7ar's. The saturation level is ±2, 
scaled by the standard deviation in 7 A p> for each event, with negative changes 
represented as red and positive changes as blue. Blues predominate, implying 
that, over time, field decay rates typically grow more positive, i.e., fields decay 
less rapidly with height. 

our CMEs. For December 12 2007 we see a slight increase in 
B,, while for May 12, 1997 we see a slight decrease in B,. A 
key point here is that CMEs can evidently occur for any type of 
evolution in B,, as with evolution in 7. Furthermore, we found 
no obvious correlation between B,{to) and either CME speed or 
flare size (see the right three columns of Table [TJ. We saw that 
PILs had a higher transverse magnetic field compared to those 



Table 3 

Evolution of Average Transverse B Field Strength at Height Closest to 42 Mm 



Initial Time" 


Final Time" 


AB, AR > b AB, AR 


< AB 


pil > AB, PlL 


<0 


-6 





5 


5 


1 


3 


-12 





6 


4 


6 


4 


-24 





6(1) 


4 


4(1) 


6 


-48 





5(1,1) 


5 


7 (2,1) 


3 


-12 


12 


6 


4 


6 


4 


-24 


24 


6(2) 


4 


4(2) 


6 


-48 


24 


5 (2,1) 


5 


7(3,1) 


3 



a Measured in hours since closest magnetogram to CME event 

Most values are not significant at a \a confidence level; the number of events for which 
values are significant are listed in parentheses with the first and second entries correspond- 
ing to lcr and 2(T significance, respectively. 



of whole ARs. This is probably due to the structure of the mag- 
netic field at the PIL: tracing field lines between two polarities, 
the potential magnetic field generally runs nearly parallel to the 
surface above the PIL, leading to a comparatively high value 
of transverse field strength. As with A7, we computed several 
AuB,'s, as shown in Table [3] There is a slight predominance 
of positive AB,pil's — five out of seven cases — which would 
have a likelihood of 0.07 if AB,pil were random, and obeyed the 
binomial distribution. As with A7AR, we calculated changes in 
normalized AB,pil for all time intervals, and plotted these dif- 
ferences in Figure [7] with the same shading and scaling. There 
is a slight predominance of blue, corresponding to a tendency 
for increases in B, . We note that there is also a predominance of 
positive AB,ar's in Table [3] but not as great as that among the 
AB,pi L 's. 

To see if increases in B, could be interpreted as due to flux emer- 
gence, we plotted the total unsigned radial flux of the entire AR 
for each model for each event in Figure[8] Error bars were com- 
puted assuming 20 Mx cm" 2 per MDI pixel, or 2.5 x 10 20 Mx 
for each ~ 12 Mm PFSS pixel. With some exceptions, there is 
a strong correlation between the evolution of the average B, at 
46.5 Mm and the total unsigned flux. While increases in B, in 
many cases appear related to flux emergence (e.g., March 13, 
2003, and Nov. 06, 2004), there are also clear counterexamples 
(e.g., Jan. 03, 2005). Also, we found examples of CMEs oc- 
curring during episodes of both increases and decreases in flux 
(e.g., Dec. 12, 2007 and Jan. 03, 2005, respectively). Observa- 
tions of CMEs while flux is both increasing and decreasing im- 
ply that changes in flux, by themselves, a re not necessary con- 
ditions for CMEs to occur. We remark that Mason & Hoeksema 
(120101) found flux emergence to be a less powerful indicator of 
large flares (cf., CMEs) than other properties of photospheric 
fields. 

We note that uncertainties for this total unsigned flux do not 
include systematic uncertain ties in the model radial ma gnetic 
fields, such as assumptions Schr iiver & De Ro"sal ((2003) made 
in converting the line-of-sight magnetic field to the radial 
magnetic field, the spatial and temporal sampling of the line- 
of-sight magnetic field measurements, Fourier ringing from the 
spherical harmonic expansion, and a well known saturation 
effect in the MDI instru ment that und erestimates flux densities 
in strong field regions dLiu et alj|200l . 



4.3. Evolution of Overlying Flux at Larger Radii 

Following Wang & Zhang (2007), we sought to investigate the 
amount of transverse flux overlying PILs. They computed 
transverse-field fluxes over low (from 1.0 to I.IRq) and high 
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Figure 6. Time evolution of the averaged transverse magnetic field B, at 46.5 Mm. The solid and dotted lines correspond to fields computed from ARs and PILs 
respectively. Some error bars and data points are slightly offset in time to avoid overlap. The solid vertical line corresponds to t = and the dashed vertical line 
corresponds to the CME time. 
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Figure 7. AB,pil between all possible pairs of B,pil's. The saturation level 
is ±2, scaled by the standard deviation in 6,pil for each event, with negative 
changes represented as red and positive changes as blue. Blues predominate, 
implying B/pil tends to increase in time. 

(from 1.1 to 1.5^0) radii crossing planes parallel to their PILs, 
as well as the low-to-high ratio of these two fluxes. We numeri- 
cally integrated the transverse field strength over radius in these 
low and high ranges for all PIL pixels, averaged the integrated 
quantities over all the PIL pixels, and took the low-high ratio. 
Each term in the ratio therefore had units of Mx cm" 1 , and can 
be thought of as the average transverse flux per pixel. 



We note that their PFSS model fields incorporated more spher- 
ical harmonic coefficients, with £ max = 225, while the default 
value of £ ma x for our models is 192. In principle, over the low 
radii from 1.0 to 1.1 ^0, higher-£ max models can more accu- 
rately represent the potential magnetic field from a given bound- 
ary condition. It is also true, however, that the corona is thought 
to be increasingly non-potential at lower heights. Consequently, 
a higher-order potential-field model does not necessarily more 
accurately represent the actual field in the low corona. To check 
if high-low flux ratios depend strongly on £ max , we computed 
the ratios with two £ max values, 192 and 225. 
In Figure [9] we plot time series of the low-to-high flux ratio, 
for two choices of the maximum spherical harmonic coefficient 
of dax, 192 (dashed) and 225 (solid). A small low-high ratio 
implies weaker low-lying transverse fields relative to transverse 
fields at greater heights. Comparing trends between events, no 
evidence for regular behavior is seen for both values of t max . 
For instance, in some events the ratio decreases (e.g., Jan. 03, 
2005), while for others it increases (e.g., June 12, 2005). 
The computations with £ max = 225 were undertaken to check 
if using higher I terms in the spherical harmonic expansion 
strongly affected the low-high ratios. We found similar evo- 
lution in the low-to-high ratio for either value of £ max : for all 
events except the May 12, 1997 event, the low-high ratios for 
each £ max are significantly correlated over our 13 sample times, 
with rank-order correlation coefficients in the range 0.4 - 0.9. 
For the May 12, 1997 event, we qualitatively compared the re- 
constructed photospheric radial fields for both values of l m3X at 
to, and found that the higher-^ radial field had more salt-and- 
pepper weak fields, which would increase the average trans- 
verse field strength at low heights. The discrepant evolution of 
the low-to-high ratio for this event does not change our conclu- 
sion that evolution in the low-to-high flux ratio does not exhibit 
any consistent trends around CME onset times. 
In addition to studying evolution of the ratio of low-to-high 
fluxes, we also considered changes in the high flux, by it- 
self. In Figure [10] we show pairwise time differences A$ of 
radially-integrated fields from 1.1 — 1.5J?©, which we denote 
$(1.1 - 1.5^0), for all possible time pairs for each event, anal- 
ogous to differences shown in Figures |5]and 7. For most events, 
these plots behave similarly to those of B,pil in Figure 7: trends 
in the magnetic field starting at 46.5 Mm continue at higher al- 
titudes up to 1.5^0. For the June 12, 2005 event, however, B, at 
46.5 Mm increases, while the integrated field from 1.1 to 1.5^0 
decreases. We note that the decreasing integrated field at high 
radii is consistent with the more-negative 7 for this event in Fig- 
ure g] 

5. DISCUSSION 

Pre-CME coronal fields are thought to be in a state of bal- 
ance between outward-directed magnetic pressure and inward- 
directed magnetic tension. Changes in the large scale magnetic 
environment of a pre-CME magnetic field could reduce the con- 
fining tension force, causing an imbalance, and thereby initiat- 
ing the CME process. Given rapid time scales of magnetohy- 
drodynamic instabilities in the corona, any significant changes 
in either outward pressure or inw ard confinement could lead 
to rapid CME onset. For instance, Kliem & Torok (2006) pro- 
posed that if the decay with radius of overlying magnetic field 
is sufficiently rapid, then the torus instability would occur. 
We looked for evidence of significant evolution in the large- 
scale magnetic environment in ten CME source regions from 
about two days before to one day after each event. We studied 
isolated active regions on the sun and the magnetic field decay 
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Figure 8. Time evolution of the total unsigned flux in each event, summed over each AR. Error bars were derived assuming uncertainties of 20 Mx cm 2 in the MDI 
magnetograms used to produce the models. The solid vertical line corresponds tat = and the dashed vertical line corresponds to the CME time. 
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Figure 9. Ratio of radially-integrated transverse field ranging between heights of 1 .0-1 . IRq to radially-integrated field between 1 . 1-1 .5Rq . Transverse field strengths 
with Lin = 225 are plotted solid while field strengths with ^ max = 192 are plotted dashed. 



rates above these regions as functions of time. We masked sub- 
sets of pixels from magnetograms into two categories: active 
regions (ARs) and polarity inversion lines (PILs). For each of 



these two categories we applied two separate methods of com- 
puting the power-law exponent, 7, characterizing the decay rate 
of transverse field with radius above the model photosphere. In 
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Figure 10. Differences of overlying integrated fields, A<E>, between all possible 
pairs of $(1.1 — 1.5Rq). The saturation level is ±2, scaled by the standard 
deviation in $(1.1 — 1.5Sq) for each event, with negative changes represented 
as red and positive changes as blue. 

Method 1, we averaged transverse magnetic field strength at a 
given radius over all pixels in the mask, then fitted the logarithm 
of these averages against the logarithm of radius to find 7. In 
Method 2, we fitted a power-law exponent 7 for the decay of B, 
with r for each pixel within a mask (AR or PIL) and then aver- 
aged these 7's. Even though the values of 7 were different for 
methods 1 & 2, a more important factor that affected the fitted 7 



was which pixels were being used in the mask: AR versus PIL. 
We then sought any relation between evolution of radial decay 
rates of B, and CME onset. On no particular time scale is there a 
statistically significant difference in the number of CME events 
whose 7 increased versus decreased. On the contrary, the ma- 
jority of 7's maintained a steady value with variations that failed 
to exceed a 1 - a confidence level. That said, while no strong 
trends were present, we did find a slight tendency for 7's to be- 
come less negative, consistent with increased transverse field 
stre n gths a t higher altitudes around CME times. 
Liu (2008) suggested transverse magnetic field strengths, at 
heights of approximately 42 Mm, may quantify the strength of 
magnetic confinement, implying that changes in this field could 
be related to CME onset. Accordingly, in addition to modeling 
7, we computed the PIL-mask- and AR-mask-averaged trans- 
verse field strengths at 46.5 Mm above our CME source regions. 
As shown in Figure [6] we also saw no common trends in any 
time scale for this magnetic parameter. 

In addition, we found that CME speeds and strengths of their 
associated flares were not significantly correlated with either 7 
or the average overlying transverse magnetic field strength at 
46.5 Mm. 

Wang & Zhang] d2007l) analyzed the ratio of transverse fluxes 
overlying flaring PILs over low (1.0-1. IT?©) and high (1.1— 
\.5Rq) altitudes. They found regions with lower ratios of low- 
to-high flux were less likely to erupt. We looked at radially in- 
tegrated transverse fields over these same ranges and found no 
evidence for significant, systematic changes in our sample, im- 
plying that no particular evolution of the low-to-high flux ratio 
is common near the times of CMEs. 

We find that evolution in the field-strength decay parameter 
7, average transverse field strength B, near 50 Mm, and radi- 
ally integrated overlying field do not follow any obvious pat- 
terns around CME onset in our potential-field models of source 
regions. Phrased another way, we find that evolution in our 
parametrizations of large-scale structure in PFSS models is not, 
by itself, a necessary condition for ejections to occur: eruptions 
were observed to occur whether our parameters increased, de- 
creased, or stayed the same. If our potential field models ac- 
curately represent the large-scale field in CME source regions, 
then our results imply that no particular large-scale evolution is, 
by itself, necessary for CMEs to occur. 

We expect, however, that large-scale field evolution could play 
an important role in the CME process, albeit one that varies 
depending upon evolution in the core fields of CME source re- 
gions. For instance, it is possible that the torus instability (or 
any other instability related to the rate of decay of field strength 
with radius) might play key roles in the CME process, but our 
results suggest that crucial evolution occurs at heights lower 
than those we modeled (primarily 46.5 Mm and above). Both 
the low spatial resolution of our PFSS models and their neglect 
of the coronal electric currents that drive CMEs, however, imply 
that we cannot sensibly extend our analysis to smaller radii. 
We remark that our models cannot accurately capture rapid 
changes in confinement. For instance, our models do not 
address the possibility of decrease in confinement via re- 
connection in the ove rlying fields, as in the breakout model 
(lAntiochos et al.|[l999l) . Investigating such processes would re- 
quire a glo bal- scale, dynamic f i eld mo del . 
Liu (2008) and lWang & Zhangld2007l) compared potential fields 
above sites of both eruptive and non-eruptive flares, and found 
that their parametrizations of large-scale confinement were 
weaker in eruptive flares. This suggested both that potential 
models could quantify confinement, and that the strength of 
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confinement is a key factor in the eruption process. We studied 
the evolution of similar parametrizations of confinement prior 
to several CMEs, and saw no clear patterns in these parame- 
ters' behavior around the time of each CME. Because we did 
not compare regions with and without ejections, our results are 
not necessarily inconsistent with previous results. We conclude, 
however, that the measures of confinement that we have inves- 
tigated are of limited utility in forecasting CMEs. 
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